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Abstract 

We investigate the non-Gaussianity of primordial curvature perturbations in the modulated 
reheating scenario where the primordial perturbation is generated due to the spatial fluctuation of 
the rate of the inflaton decay to radiation. We use the 5N formalism to evaluate the trispectrum 
of the curvature perturbation as well as its bispectrum. We give expressions for three non-linear 
parameters /jvx, tjvl and g^L hi the modulated reheating scenario. If both the intrinsic non- 
Gaussianity of scalar field fluctuations and third the derivative of the decay rate with respect to 
the scalar fields are negligibly small, g^L has at least the same order of magnitude as /jvl- We 
also give a general inequality between /tvl and tjvz,, which is true for other inflationary scenarios 
as long as the primordial non-Gaussianity comes from super-horizon evolution. 

PACS numbers: 98.80.Cq 
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I. INTRODUCTION 



Recent observations of cosmic microwave background (CMB) anisotropies give strong 
evidence that primordial density fluctuations are almost Gaussian, scale-invariant, and adi- 
abatic [lj]. During inflation, they are generated as vacuum fluctuations of light fields and 
are stretched to cosmological scales to explain the large scale structure of the universe 0] • 
Such a light field responsible for density fluctuations has been considered to be the inflaton 
itself for a long time. 

Recently, alternative candidates for such a light field have been proposed. One attractive 
example is the curvaton [l, 0] , which is effectively massless and acquires fluctuations during 
inflation. After inflation, it becomes effectively massive and contributes to a non-negligible 
fraction of the energy density of the universe. Then, after it decays, eventually density 
fluctuations induced by the curvaton are converted to adiabatic ones and can dominate 
over those generated by the inflaton itself. Another interesting candidate is a light field 
whose expectation value determines the coupling constant of the inflaton to standard model 
particles j5[. Such a light field will fluctuate during inflation, which leads to fluctuation 
of the decay rate of the inflaton. Then, this latter spatial fluctuation induces that of the 
reheating temperature, which eventually generates curvature perturbations. Some variants 
of these two alternative candidates have also been considered related to the fluctuations of 
masses and annihilation cross sections and the preheating mechanism [tJ. 

In order to determine which light field is actually responsible for primordial density fluc- 
tuations, the deviation from Gaussianity of the curvature perturbations is of great use. If 
density fluctuations are completely Gaussian, their bispectra, characterized by a parameter 
Jnl, and (connected) trispectra, characterized by parameters g^L an d tjvl vanish. There- 
fore, estimation of the bispectra and trispectra is important to identify the light field. In Q|, 
it was shown that the bispectrum is significantly suppressed by the slow-roll parameters up 
to an undetectable level in a single field slow-roll inflation model. 1 After 0], the trispectrum, 
as well as bispectrum, were evaluated in a multi-field configuration and a curvaton scenario. 
Though spectra are still suppressed by slow-roll parameters in a multi-field configuration 
0,[n|, they can be significantly large in the curvaton scenario [12J, Il3f . On the contrary, 



only the bispectrum is calculated in the modulated reheating scenario [14|, [15j, [16J. Since the 
bispectrum can be large in both the curvaton scenario and the modulated reheating scenario, 
trispectrum may be useful to discriminate between them. Though the present constraint on 



the trispectrum is not particularly severe, and is roughly given by \tnl\ < 10 8 17] . a value 
of \t~nl\ ~ 560 will be detectable by the Planck satellite |l8|]. 2 

The main purpose of our paper is to estimate the trispectrum in the modulated re- 
heating scenario. We use the 5N formalism, which is a powerful approach to evaluate the 
non- Gaussianity of curvature perturbations simply because it requires a homogeneous back- 
ground solution @, Hi], [H, H3] - Then, we use the same formalism to evaluate the trispectrum 



of curvature perturbations as well as its bispectrum. 

This paper is organized as follows. In the next section, we give a brief review of the SN 



1 It has recently pointed out that there is a possibility that an all-sky 21-cm experiment is sensitive to a 
value of \f NL \ ~ 0.01 fl. 

2 In order to obtain these constraints, one must compute not only the non-Gaussianities of primordial 
fluctuations but also other non-linear effects entering in the CMB such the non-linear evolution of the 



perturbations after inflation 19 1 
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formalism and a definition of the three non-linear parameters Jnl, t~nl and g^L given in 



13l . 1 171 l24j . We also give a general inequality between f^L and tnl which is not found in 
the literature. In Sec. Ill, we study the background dynamics in the modulated reheating 
scenario. In Sec. IV, we study the perturbations in this scenario and give expressions for 
Jnl, tnl and g^L- The final section is devoted to a summary. We use the units 811G = 1. 



II. 5N FORMALISM 



According to the 5N formalism [20|, |2l|, |22|, |23j, the curvature perturbation on a uniform 



energy density hypersurface ( at time tf is, on sufficiently large scales, equal to the pertur- 
bation in the time integral of the local expansion from an initial flat hypersurface (t = ti) 
to the final uniform energy density hypersurface. On sufficiently large scales, the local ex- 
pansion can be approximated quite well by the expansion of the unperturbed Friedmann 
universe. Hence 

C(tf,x) = N(ti,tf,x) — (spatial average), (1) 

where the e-folding number N(ti,tf,x) is defined by the time integral of the local Hubble 
parameter, 

ftf 

N(ti,t f ,x) = / H(t,x)dt. (2) 
Ju 

In many inflationary scenarios, which include the modulated reheating scenario, the dy- 
namics of the universe between ti and tf is determined by the values of the relevant scalar 
fields (p 1 at ti and by the e-folding number which becomes a function of 7 (t i; x). Hence the 
curvature perturbation at tf is given by 

((t f , x) « N I S(f) 1 + ^A/j50 7 50 J H (spatial average), (3) 

where 5(f) 1 is the perturbation of (j) 1 on the flat hypersurface at ti, and Nj, Njj, ■ ■ ■ are given 
by 

_ j dN AT d 2 N 

Ni = w Nij =ww> ••- (4) 

Because solutions of the unperturbed Friedmann equation give Ni, Njj, ■ ■ -, the knowledge 
of the background solutions is enough to know the higher order correlation functions of (. 
The connected parts of the power spectrum, bispectrum and trispectrum are defined as 

(C^ k2 ) c ^(27rfP < (k 1 )S(k 1 + k 2 ), (5) 
(CfeCfeCfi,)* = (^) 3 B ( (h, k 2 , k 3 )6(h + k 2 + fc s ), (6) 
(CfcCftCftCs^e = (2vr) 3 T c (^, k 2 , h, h)S(h + k 2 + h + h). (7) 

Here (■ ■ -) c means that we take the connected part of (• • ■). Using Eq. ([3]), we can express 
P{, and in terms of the correlation functions of 5(f) 1 , which are given to leading order 
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by m 

P c (k)=NjNjP IJ (k), (8) 
B c (h, k 2 , k 3 ) = NiNjNkB ijk (k± , k 2 , k 3 ) + NjNjkNk (P ik \h)P JL (k 2 ) + 2 perms) (9) 
T c (h, fc 2 , fc 3 , k 4 ) = NjNjN K N L T IJKL (k u fc 2 , fc 3j fc 4 ) 

+N U N K N L N M (P IK (h)B JLM (k 12 , k 3 , h) + 11 perms.) 

+iV /J 7V Ki iV M 7V JV (P Ji (fc 13 )P /M (fc 3 )P J7V (A; 4 ) + 11 perms.) 

+N IJK N L N M N N (P IL (k 2 )P JM (h)P KN (h) + 3 perms.) , (10) 

where kij = \ki — kj\ and P IJ , B IJK , T IJKL are the power spectrum, bispectrum and 
trispectrum of the scalar fields respectively, defined by 

= (2 7 r) 3 P IJ (k 1 )5(k 1 + fc 2 ), (11) 

£^ 2 ^|) c = (2tt) 3 j B w (A ;i , fc 2 , fc 3 )5(£ a + £ 2 + fc 3 ), (12) 

£^ 2 ^|^ 4 ) c = (27r) 3 T /Ji ^ i (A ;i , fc 2 , fc 3 , fc 4 )5(£i + fc 2 + fc 3 + k). (13) 

If 5c/) 7 are independent Gaussian variables with the same variance which we denote as P, then 
P IJ = PS IJ and B IJK and T IJKL vanish. In such deviation from Gaussianity of the 

primordial perturbation comes only from super-horizon evolution, and the non-Gaussianity 
is characterized by three constant parameters /ml, tnl and qnl defined by 

B c (h, k 2 , k 3 ) = p NL (P^k^P^h) + 2 perms.) , (14) 

T f (fci, k 2 , k 3 , k 4 ) = r NL (Pc(A;i 3 )P c (A; 3 )P c (A; 4 ) + 11 perms.) 

54 

+^9nl (P<(k 2 )P c (k 3 )P c (k 4 ) + 3 perms.) . (15) 

Using Eqs. (1H1)-( |T0|) . we have the following expressions, 

5 NjNjN ij 

fNL ~ e^v^f' (16) 

NjjN ik N j N k 
TNL = (N L N^ ' (17) 



#iVL = — — „ 3 • (18) 



25 NhkN^'N* 
54 (N L N L f 

Eg . f[16l) has been given in 24j. Eo. (fTTl) has been given in the arXiv version of 17] and also 
in [l3j. Eq. ffl8|) has been given in [13j. 

Here we provide a relation between /jvl and t/vl which is not found in the literature. 
From the Cauchy-Schwarz inequality, we have the following 

36 „ 9 . 
t^l > -f NL . (19) 

We have equality if and only if the vector Ni is an eigenvector of the matrix Njj. The single 
inflation model yields tnl = ff/ivx- However in multi-field inflation, there is a possibility 
that two vectors Nj and NijN j are nearly orthogonal. In such a case, /jvl is very small 
but tjvz, and possibly also g^i remains finite. Hence the leading non-Gaussianity comes not 
from the bispectrum but from the trispectrum. 
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III. BACKGROUND DYNAMICS OF THE MODULATED REHEATING SCE- 
NARIO 



In the modulated reheating scenario 0, the decay rate T of the inflaton S is a function 
of scalar fields (j) 1 (not necessarily a single field) which are light during inflation. We assume 
that fluctuation of the inflaton field generates negligible curvature perturbation. Then the 
detailed form of the inflaton potential U(S) during inflation is not important for the scenario. 
We only require that U (S) around the minimum is approximated well by a term quadratic in 
S. After inflation, the inflaton oscillates around the minimum of the potential. The energy 
density of the inflaton ps averaged over one period of the oscillation behaves as a function 
of the e-folding number oc e~ 3N . Hence we regard ps as dust. The inflaton decays into 
radiation with rate T which depends on the expectation values of </> 7 . Then the background 
equations are given by 

^ + 3p s = -I„ s , (20) 

H 2 = i(p s + Pr), (22) 

where p r is the energy density of radiation. Spatial fluctuations of these light fields induce 
the fluctuation of the decay rate and the curvature perturbation. By solving the above 
equations from the end of inflation to the completion of reheating with the initial conditions 
Ps(0) = po = 3i?o, Pr(0) = 0, we can obtain a relation between N and T. Until the Hubble 
function drops to Hf (Hf <C T), the e-folding number N can be written formally as 

N=ho e ^ +Q (l), (23) 



Hf \H 



where 

exp 



oc 



AQ(T/H ) ee / dN> Jf^/"'^ 1 - (24) 



o 



For the two limiting cases, we have the approximate form of Q(x) as 

Q( x ) = — + 0{x~ 2 ) x > 1, (25) 

Q(x) = -^\ogx + 0(x) a;<l. (26) 

For arbitrary x, we do not have an analytic form for Q(x) and we have to solve the back- 
ground equations numerically. We show Q(x) calculated numerically in Fig. [TJ We also find 
an accurate fitting formula for Q(x) of the form, 

W)=lfe(l + -). (27) 
4 r(x) + 3 \ x J 

For r(x) = 2.16x ' 72 , the relative error is within 2 percent. For r(x) = 1.7x ' 9 + 0.3x ' 18 , the 
relative error is within 0.5 percent. 
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FIG. 1: Plot of Q(x) calculated numerically. The fitting formula Qat(x) lies within the solid line 

IV. PERTURBATION IN THE MODULATED REHEATING SCENARIO 

From Eq. fl23l) . we can calculate 5N in terms of 5T, 

6N = xQ'(x) — + -Q"(x) ^— J + -Q"'(x) ( — 



+ 



(28) 



In the modulated reheating scenario, the decay rate depends on the scalar fields which are 
almost massless during inflation. Hence the perturbation of the decay rate can be written 
as a function of the perturbation of the scalar fields, 



5T = r 7 c 



IJKt 



2 T T 6 
From these equations, Nj, N TJ , N IJK can be written as 



+ 



N I = xQ'(x)^., 



N; 



UK 



xQ'(x 



r + ^ 

^IJK 

r 



r, 






r 


r ' 




(x) 


( T r 


r 'jk 




r 



(29) 

(30) 
(31) 



+ — — + 

r r 



r K Tjj 



+ x*Q">(x)^ L -p2) 



Substituting these latter into Eqs. (l8|)-(fT0l) yields the power spectrum, bispectrum and 
trispectrum of the primordial curvature perturbation in the modulated reheating scenario. 
Note that Eqs. (130]) - (132]) are only correct if tj is taken as the time when inflation ends. Hence 
in Eqs. (l8l- (fTUl) . we must use the power spectrum, bispectrum and trispectrum of scalar field 
fluctuation at the end of inflation. 



Zaldarriaga evaluated the bispectrum B of the curvature perturbations at reheating 
by taking into account not only the nonlinear evolution of the scalar fluctuations but also 
their intrinsic non-Gaussianities [13]. In particular, in order to estimate the intrinsic non- 
Gaussianity of a scalar field, he calculates the leading order quantum three-point correlation 
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function with the assumption that T depends only on a single field 0, 



(<PkM^ 2 (v)^ 3 (v))c = -i(0\ / drf [^ 1 (r?)^ 2 (r ? )^ 3 (r ? ),^(r ? ')]0), (33) 



where 77 = J_ oc — is the conformal time. In Ref. 14], mode functions in pure de-Sitter 
spacetime are used throughout the evolution of the above correlation function, though the 
use of such mode functions are not necessarily justified during the oscillatory phase after 
inflation. In our formalism, we have only to know the field fluctuations at the end of 
inflation. However, even in this case, mode functions in pure de-Sitter spacetime are not 
good approximations after the relevant scale crosses the horizon. This is because the Hubble 
parameter may decrease by an order of magnitude from the horizon exit to the end of 
inflation. In this paper, instead of using Eq. (133]) until the end of inflation, we use it 
only until the relevant mode exits the horizon. After the horizon crossing, we evolve 5(f) 1 
by perturbing the classical unperturbed equations, in keeping with the spirit of the 5N 
formalism, which correctly takes into account the evolution of the Hubble parameter. 

We start from the evolution of scalar field fluctuations on super-horizon scales, which can 
be described well by the classical treatment. Then, the background equations are given by 

dV ( 1 dH\ dS 1 V 1 



dN> + y 3+ HjN)jN + lP =0 ' (34) 

We assume that scalar fields slow- roll during the whole epochs of interest, which enables us 
to approximate the background equations as 



V 1 



dN~ 3 + f# 2 



(35) 



Throughout this paper, it is also assumed that the total energy density of the universe is 
dominated by the inflaton potential, and the dependence of cf) 1 on H is negligible. 

Let JV* be the e-folding slightly after the horizon crossing and Nf be the e-folding at the 
end of inflation. Slightly after the horizon crossing, the scalar field fluctuations turn into 
classical variables with their magnitude given by 5(f>l. Then 5(f) 1 (N) after N* is a function 
of 5(fil, which can be Taylor-expanded with respect to 5(f>{. For the sake of the evaluation of 
the leading bispectrum and trispectrum of scalar fields at N = Nf, it is enough to expand 
5(f) 1 (N) to third order in 5<p^. Up to third order, Eq. (135!) yields 

tf{N f ) = A r j(N f , N.Wi + \& JK {N h N m )5rf5tf + \^ JKL {Nf, NjStfStffStf, (36) 
where 

Aj(N,N') = [Texp ( [ dN" P{N")) Y , (37) 



A" 



r N f 

& JK {N f , AQ = / dJVA£(JV), N>)Q L MN (N')Af(N', N*)A%(N>, N*), (38) 



% KL {N f , N.) = - dN' A^N, N^QiMN^AfiN', N.) 
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x / dN" A»,(N', N")Q N PQ (N")A P K {N\ N,)Af(N\ NJ 

T> A I (AT AT'Wl' fAT>\\J'fAT> AT \ A K> ( A T> AT \ A L' ( AT' 



+ / dN' Aj,(N, N )Rj, k , l ,{N')Kj (N 1 , N.)A% (N', N*)A h L (N', N*p9) 



i v! 



Pj = -7TTWlk> ( 4 °) 



3 + f # 2 
1 V 



i 



= "JTf #■ (41) 

1 V 7 

RjKL = - 3+ H> jj 2 L - (42) 

A. Power spectrum of 5c/) 1 at the end of inflation 

To leading order, the two-point function of 5(f) 1 at the end of inflation is given by 

P IJ (N f ) = A T K (N f , N*)A J L (N f , N*)P? L (h). (43) 
Here, slightly after the horizon crossing, P/ J is given to leading order by, 

yu _ ( 27 3 ( H *\\u _ D xu 



x-*&W JsP >*- (44) 



Any corrections to Eq. (jUJ) are suppressed by the slow-roll parameters [13j. Only fluctuations 
that extend to super-horizon scales at the end inflation contribute to the fluctuation of V. 
Such fields must be massless during inflation and we assume A 1 « 5 1 , which yields the 
corresponding power spectrum as 



P u {N f )^PJ u . (45) 
B. Bispectrum of 5(f) 1 at the end of inflation 

To leading order, the bispectrum of 5(f) 1 at the end of inflation is given by 
B IJK (h, k 2 , k 3 )(N f ) ~ Bi JK (k u k 2 , k 3 ) + Q IJK (N f , TV*) (P^hjP.fa) + 2 perms.) . (46) 
E>l JK , which is the bispectrum slightly after the horizon crossing, can be evaluated by 



quantum perturbation theory [14 



Bi JK (h, k 2 , h) = -2V IJK Re ig kl (v*)9k 2 (v*)9k 3 (v*) I dij a 4 (rj)g* kl {rj)g* k2 (rj)g* k3 (r]) , (47) 



4/ 



oo 



where gk{ff) is the mode function in de-Sitter spacetime and is given by, 

g k (v) = ^£^ 2 (l+ikr,)e- ikr >. (48) 
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Soon after the horizon crossing, kit] becomes smaller than unity. In such a phase, the leading 
bispectrum of Eq. (|47|) is given by 



Bi JK (k u k 2 , h) = Vj 



El k\ 

3 V UK 



4 (hk 2 k 3 

7 + I log M A (-1 + - k ^P) + t- %^1(,49) 



2 taV t ; * y y V 3 fcf k* J ' 9 k 2 t 



where 7 = 0.577... is Euler's constant. 3 Because the second line in Eq. (1491) is 0(1), 
= 0{Vuk/ Hf). On the other hand, the bispectrum of the second term in Eq. ( )46|) can 
be estimated as 

e IJK (N f , njp.mpm = o (Jp £ f dN ^wy) ' (50) 

which is expected to be larger than 5* because Nf — N* is typically O(10). Hence the 
leading bispectrum of <50 7 at the end of inflation comes from super-horizon evolution, 

B IJK {N f ) « e IJK (N f , N*) (P^kjP^h) + 2 perms.) . (51) 

Thus, the corresponding non-linear parameter f^L becomes 4 

e f 1 rr f e' 1 nyffe q" (x) 

5 JNL xQ'(x)T K T K xQ'(x)(T K T K f 2 Q' 2 (xY 1 1 

where Q 1 and FL) are projected vectors of Q IJK and r /J , respectively, 



e 1 = y^Q ijk , (53) 

(2) {VKV Kf/2- 



r? 2 ) = 7-4rr^- (54) 



One should notice that the first term in Eq. (1521 represents intrinsic non-Gaussianity due to 
cubic interactions of the scalar fields. The second term comes from non-linearity between V 
and cj) 1 . The third term comes from non-linearity between ( and 5T, which only depends on 
T through the argument x = T/Hq of Q(x). Q"(x)/Q' 2 (x) takes minimum value 6 at x = 
and becomes larger for x > 1 (see Fig. [2]). Hence if both the intrinsic non-Gaussianity of the 
scalar field fluctuations and the non-linearity between T and 7 are negligibly small, then 

InL > 5. 



3 This expression for the bispectrum is slightly different from that given in Rcf. 14j, though the difference 
is not essential. Thanks to a second look by the author of Ref. 14j, we reach the same result, which is 
given above. 



Note that the definition of /jvl here is different in sign from that in 16 1 
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10° 



10" 



-itr 



10' 



10' 



10" 




FIG. 2: Plot of Q"(x)/Q ,2 (x) 



C. Trispectrum of 5c/) 1 at the end of inflation 

To leading order, the trispectrum of 5(f) 1 at the end of inflation is given by 

T IJKL {k u k 2 , k 3 , h)(N f ) = T! JKL (h, k 2 , k 3 , h) 

+ef M (Nf,N*)ej M (N f ,N*) (P*(h)P4h)P4k 13 ) + ll perms.; 

perms. J , 



(55) 



As in the case of the bispectrum, the leading contribution to Tl JKL can be evaluated as [l 

T! JKL (h,k 2 ,k 3 ,h) = -: 



'IJKL 



xRe 

For kiT]* <C 1, Eq. fl56l) reduces to 
T! JKL (h,k 2 ,k 3 ,k 4 



.(56) 



H 4 kf 



8(kik 2 k 3 k4 



-3V1JKL 



x 



-log(/^*) 2 + 7 



1 

3 



kf 



4 Y.i<jkikj k x k 2 k 3 k± 
+ 
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kf 



(57) 



where k t = k\ + k 2 + k 3 + k^. T* is typically smaller than the other terms on the right 
hand side of Eq. ( |55i) . which corresponds to the super- horizon evolution of the trispectrum 
due to the fourth order interactions. Hence, as in the case of the bispectrum, the leading 
trispectrum of 5(f) 1 at the end of inflation comes from super-horizon evolution, 



T IJKL (h,k 2 , k 3 , h)(N f ) = Qf M (Nf, N*)Qj M (Nf, N*) (P*(fci)P*(A; 2 )P*(fc 13 ) + 11 perms.) 



+Zf JK (N f , i\y (P^P^P^h) + 3 perms.) 



(58) 
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FIG. 3: Plot of Q'"{x)/Q ri {x 



Therefore, the corresponding non-linear parameters tnl and g^L are given by 

i r 2 e J ej 2r 2 f( 2) e / 2x 2 Q"(x)rr / e / 



tnl 



(xQ'(x)) 2 T K T« ( x Q'( x )) 2 (T K T^f 2 (xQ>(x)) 3 T K TK 



r 2 f( 2) f /(2) 2x 2 Q"(x) rr 7 f( 2) g«, 



x 



L\2 



54 

25 9NL 



(xQ'(x)Y(T L T L ) 
T 2 E 



+ 



[xQ'{x)f (r K r«f 2 

(2) 



Q' 4 (x) ' 

r 2 fLe 7 x 2 Q"(x) rr / e / 



+ 



(xQ'ix^TjT 1 (xQ'(x)) 2 (TjT J f 2 (xQ'(x)) 3 T K T K 

i r 2 r 7 ff 3) 3x 2 Q"(x) rf( 2) r 7 g- (x) 



where 



(xg'(x)) 2 (r^) 2 



TiTjT k T l IJKL 



(xQ'(x)) 3 (TjTJf 2 Q ,3 (x) 



(T M T M ) 2 



r 1 

1 (3) 



TjT K r IJK 

TrT L 



(59) 



(60) 



(61) 



As in the case of the bispectrum, the last term of g^L represents the contribution from the 
non-linear relation between £ and 5T and depends on T only through the argument x of 
Q(x). Q'" (x) j Q' 3 (x) takes minimum value 72 at x = 0, and larger values for x > 1 (see 
Fig. [3]). We see that all three non-linear parameters, /jvl, ttvl and <7jvx, can be written only 
in terms of the vector quantities T 1 , T^, and 9 7 . 

If r JJX and the intrinsic non-Gaussianity of scalar field fluctuations are negligibly small, 
we have the following relation 



5 Q"(x) 
9NL -3Q^) fNL ~ 

In particular, for x <C 1, this relation reduces to 



25 Q"'(x) 
54Q' 3 (x) 



9nl = 10/. 



NL 



50 



(62) 



(63) 
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Hence g^L has the same order of magnitude as Jnl for i < 1 and becomes much larger than 
Jnl for x ^> 1. This is in contrast with the case of standard single-field inflation, where g^L 
is to second order in the slow-roll parameters. Notice that the relation Eq. (1621) depends on 
T only through the argument x = T/H and is independent of the functional form of T(0) 
except for conditions on T IJK . 



V. CONCLUSIONS 

The modulated reheating scenario generates a primordial curvature perturbation due 
to the spatial fluctuations of the inflaton decay rate. Such a scenario induces larger non- 
Gaussianity in the perturbation than that of the simple inflation scenario. Future observa- 



tions such as the Planck satellite [18] may be able to detect or constrain the level of the 
non-Gaussianity, and thereby distinguish the different scenarios. 

We have given expressions for the power spectrum, bispectrum and trispectrum to lead- 
ing order in the modulated reheating scenario, allowing for multi-field dependence on the 
inflaton decay rate. The leading contribution to the bispectrum and trispectrum comes from 
the super-horizon evolution. If the intrinsic non-Gaussianity of the scalar fields and third 
derivative of the decay rate T IJK are subdominant, then we have a simple relation between 
two non-linear parameters /jvl and gNL, which is independent of the detailed form of the 
decay rate except for the condition on T IJK . Hence, g^i has at least the same order of 
magnitude as /tvl- 

We have also given a general inequality between the bispectrum and the trispectrum 
t~nl > ff/ivL which is true for other inflationary scenarios as long as the non-Gaussianity 
comes from the super- horizon evolution. This inequality allows J'nl to be vanishingly small 
while tnl remains finite. In such a case, the leading non-Gaussianity comes not from the 
bispectrum but from the trispectrum. 
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